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The problem of anomalous diffusion in momentum space is considered for 
plasma-like systems on the basis of a new collision integral, which is appropriate 
for consideration of the probability transition function (PTF) with long tails in 
momentum space. The generalized Fokker-Planck equation for description of 
diffusion (in momentum space) of particles (ions, grains etc.) in a stochastic system 
of light particles (electrons, or electrons and ions, respectively) is applied to the 
evolution of the momentum particle distribution in a plasma. In a plasma the 
developed approach is also applicable to the diffusion of particles with an arbitrary 
mass relation, due to the small characteristic momentum transfer. The cases of 
an exponentially decreasing in momentum space (including the Boltzmann-like) 
kernel in the PT-function, as well as the more general kernels, which create the 
anomalous diffusion in velocity space due to the long tail in the PT-function, are con- 
sidered. Effective friction and diffusion coefficients for plasma-like systems are found. 

PACS number(s): 52.27.Lw, 52.20.Hv, 05.40.-a, 05.40.Fb 

I. INTRODUCTION 

Diffusion in coordinate and in momentum (velocity) space is of fundamental importance 
and has attracted a growing interest during many years, since the description of this processes 
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provides a simplified and effective key for addressing of many problems of kinetic theory 

Deviations from a linear time-dependence, i.e. < r 2 (t) >~ t, of the mean square displace- 
ment in coordinate space has been experimentally observed, in particular, under essentially 
non-equilibrium conditions or for some disordered systems. The average square separation of 
a pair of particles passively moving in a turbulent flow grows, according to Richardson's law, 
with the third power of time [1] . For diffusion typical for glasses and related complex systems 
[2], the observed time dependence is slower than linear. These two types of anomalous dif- 
fusion are obviously characterized as superdiffusion < r 2 {t) >~ t a (a > 1) and subdiffusion 
(a < 1) [3]. For a description of these two diffusion regimes a number of effective models 
and methods have been suggested. The continuous time random walk (CTRW) model of 
Scher and Montroll [4], leading to subdiffusion behavior, provides a basis for understand- 
ing photoconductivity in strongly disordered and glassy semiconductors. The Levy-flight 
model [5], leading to superdiffusion, describes various phenomena as self-diffusion in micelle 
systems [6], reaction and transport in polymer systems [7], and is applicable even to the 
stochastic description of financial market indices [8]. Both cases can be effectively described 
by generalized diffusion equations with fractional derivatives in time or in coordinates, re- 
spectively [9] . For example, different aspects of the anomalous diffusion in coordinate space 
were considered within this scheme in Refs. [10], [11]. 

However, recently a more general approach has been suggested in [12], [13]. This one per- 
mits to reproduce the results of the standard fractional differentiation method in coordinate 
space, when the latter is applicable, and enable to describe more complicated cases of anoma- 
lous diffusion processes. In [14] this approach has been applied also to the case of diffusion 
in a time-dependent external field in coordinate space. In what follows, we concentrate on 
the problem of diffusion in momentum space in application to plasma systems. 

Problems of diffusion in momentum space have been considered for plasmas in the fun- 
damental study by Landau [15] and later by Rostoker and Rosenbluth, Lenard, Balescu, 
Klimontovich and many others. Various theoretical and experimental aspects of these in- 
vestigations can be found in [16-18]. 

Here our main interest is focused on anomalous diffusion in momentum space by using 
the methods developed in [12], [13] for coordinate space. Recently, see [19], a new kinetic 
equation for anomalous diffusion in velocity space has been derived on the basis of an 
appropriate expansion of PTF and some particular problems were investigated on this basis. 
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In the present paper the problem of anomalous diffusion in momentum (velocity) space will 
be considered for plasma-like systems. 

Some aspects of anomalous diffusion in velocity space have been investigated during 
the last decade in a number of studies [20-23]. In particular in [22] the phenomenological 
equation for anomalous diffusion in velocity space for a magnetized plasma has been obtained 
on the basis of the Langevin model with the linear friction proportional to the particle 
velocity and with non-Gauss noise. The corresponding Fokker-Planck type equation included 
a diffusion term with a fractional derivative and the usual drift term with the first momentum 
derivative. The applicability of this equation is limited, due to phenomenological nature of 
the considered model. We should mention the analogy of this equation with the structure 
of the equation for anomalous diffusion in coordinate space under action of an external field 
[11], [14]. An equation similar to that in [22] has recently been applied in [23] to describe the 
evolution of the velocity distribution function of strongly non-equilibrium hot and rarified 
plasmas. Such types of plasma exist, e.g., in tokomaks. 

On the whole, in comparison with anomalous diffusion in coordinate space, anomalous 
diffusion in velocity space is has been studied to a modest extent. 

In this paper the problem of anomalous diffusion in momentum (velocity) space will be 
considered for plasma-like systems. In spite of formal similarity, the physical (and math- 
ematical) nature of diffusion in momentum space is very different from that in coordinate 
space. This is clear already from the fact that momentum conservation, which takes place 
in momentum space, has no analogy in coordinate space. 

Diffusion in velocity space for the cases of normal and anomalous behavior of the PT- 
function is presented in Section II. Starting from the argumentation based on the Boltzmann 
type of the PTF, we describe the new approach to the kinetic equation, which in fact can 
be applied to the wide class of PTF functions based on the prescribed distribution function 
for one (light) sort of particles. The anomalous diffusion in momentum space for plasma 
is analyzed in Section III on the basis of the Boltzmann-type kernel for PTF. Models of 
anomalous diffusion for plasma-like systems are considered in Section IV. In Section V the 
generalized Fokker-Planck equation for diffusion, written for the formal Fourier-component 
f(s,t) of the distribution function f(p,t), is represented in partial derivatives in velocity 
space. This representation is possible only in particular cases of the power dependence of 
the coefficients in the generalized diffusion equation. 
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II. DIFFUSION IN THE VELOCITY SPACE ON THE BASIS OF A 

MASTER- TYPE EQUATION 

Let us consider now the main problem formulated in the introduction, namely, diffusion 
in momentum space (V-space) on the basis of the master equation, which describes the 
balance of grains entering and leaving point p at time t (see, e.g. [24,25]) 

df(p,t) 



dt 



J dq{w(p,p')f(p',t)-w(p',p)f(p,t)}. (1) 



The structure of this equation is formally similar to the master equation (see, e.g. [13]) in 
coordinate space. Here w(p, p') is the kernel describing the transition probabilities. Note 
that there is only one rather general condition which w(p, p') should satisfy if the stationary 
solution exists: the balance condition of the detailed for a stationary distribution function 
f st (p), which reads 

w(P.P') = / 8f (p) (2) 
w(p\p) f st (p') { ] 

In the following analysis we use a form of the master equation [26] equivalent to ([[]) 

df(p,t) 



dt 



J dq{W(q,p + q)f(p + q,t)-W(q,p)f(p,t)}. (3) 



The probability transition H^p, p') describes the probability for a grain with momentum p' 
at point p' in momentum space to transfer from this point p' to the point p per unit time. 
The momentum transfer is equal q = p' — p. Of course, as mentioned above, the overall 
momentum has to be conserved. 

Assuming in the beginning that the characteristic changes in momentum are small one 
may expand Eq. ([3]) and arrive at the Fokker-Planck form of the equation for the density 
distribution f(p,t) 



df(p,t) d 



A a (p)f{p,t) + -^-{B a p(p)f(p,t)) 



(4) 



dt dp a 

A a (p) = J d r qq a W(q,p); B a p(p) = i J d r qq a q p W{^ p). (5) 

The coefficients A a and B a p describe the friction force and diffusion, respectively. Here r is 
the momentum space dimension. 

Because the velocity of heavy particles is small, the p-dependence of the PTF can be 
neglected for the calculation of diffusion coefficient, which in this case is constant B a p = 
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SapB, where B is the integral 

B = ±- / d r qq 2 W{q). (6) 



2r 

If we neglect the p-dependence of the PTF at all, we arrive to the coefficient A a = 
(while the diffusion coefficient is constant). In this approach, which is known to be incorrect, 
the coefficient A a for the Fokker-Planck equation can be determined on the basis of the 
argument that the stationary distribution function is Maxwellian. In this way we arrive at 
the standard form of the coefficient MTA a (p) = p a B, which is one of the forms of Einstein's 
relation. For systems far from equilibrium this argument is not acceptable. 

Following [1], we now generalize, the Fokker-Planck approach to find the coefficients of the 
kinetic equation, which are applicable also to slowly decreasing PT functions. We apply a 
more general approach, based on the difference of the velocities for light and heavy particles. 
For calculation of the function A a we have take into account that the function W (q, p) is 
scalar and depends on the variables q, q • p,p. Expanding W(q, p) on q • p one arrives at 
the following approximate representation of the function W(q, p): 

W(q, p) ~ W(q) + W'(q)(q ■ p) + \w"{q){q ■ p) 2 , (7) 

where W'(q) = dW{q, q • p)/<9(qp) | q . p=0 and W"{q) = d 2 W(q, q • p)/<9(qp) 2 | q . p=0 . 
Then, with the necessary accuracy, A a equals 

A a (p) = J d r qq a q m W'(q)=p a J d r qq a q a W\q) = ^ j d r qq 2 W'(q) (8) 

If the function W(q,p) satisfies the equality W\q) = W(q)/2MT, we arrive at the usual 
Einstein relation 

MTA a (p)= Pa B (9) 

Let us check this relation for Boltzmann-type collisions, which are described by the PT- 
function W(q,p) = w_e(q, p) [13]: 



2ir [°° _ da 



q 2 

arccos (1 - - ), u 



f b (u 2 + v 2 -q-v/fi), (10) 



where (p = Mv) and da(x, u )/d° , A* an d fb are the differential scattering cross-section, 
the mass and distribution function for the light particles, respectively. In Eq. ( jTOi) we 
took into account the approximate equalities for the scattering of light and heavy particles 
q 2 = (Ap) 2 = p' 2 (l — cos6) and 9 ~ \i where p' = fiu is the momentum of the light particle 
before collision. 
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For the equilibrium Maxwellian distribution /° the equality W'(q) = W(q)/2MT is 
evident and we arrive at the usual Fokker-Planck equation in velocity space with constant 
diffusion and friction coefficients D = B /M 2 and friction /3 = B/MT = DM/T, respectively, 
which satisfy the Einstein relation. 

For some non-equilibrium situations the PTF, as a function of the variable q, possesses a 
long tail. In this case we have to derive a generalization of the kinetic equation in the spirit 
of the analysis of the coordinate case [13,14], because the diffusion and friction coefficients in 
the form of Eqs. ([6]) and (JHJ) diverge for large q if the functions have an asymptotic behavior 
W(q) ~ l/q a with a < r + 2 and (or) W'(q) ~ 1/q? with /3 < r + 2. 

Inserting the expansion (J7|) for W(q, p) in Eq. ([3]) (as an example we choose r = 3; the 
analysis for arbitrary r runs in a similar way) we arrive at a new collision term of the kinetic 
equation, which can be considered as a generalization of the Fokker-Planck equation for 
anomalous diffusion in velocity space [19]: 



4fM ^A{s)f{s) + B a {s) dfM 



(11) 



dt ds a 
In fact, as shown in [19], in the expansion ([?]) for W(q, p) we have to keep (with the 
necessary accuracy) only the terms linear in qp and p. The function /(s) in ( TTTT) is the 
Fourier-component /(s) = J exp(zps)/ 9 (p, t) and the coefficients are equal to 



A(s) 



dq[exp(-z(qs)) - 1]W(?) = 47r 



poo 

\ dqq 2 


sin (qs) 


to 


qs 



W(q), 



(12) 



B r . 



8 n B(s): B(s) 



s 

4tt roc 



dqqs[exp(— z(qs) — 

sin(qs) 



dqq 



o 



cos (qs) 



qs 



W'(q). 



(13) 



Here we took into account the existence of the small parameter yu/M and we omitted the 
small on this parameter terms of order p 2 and W" in Eq. ([7]) . 

For the isotropic function f(s) = f(s) one can rewrite Eq. (ITTj) in the form 

(14) 



df ^ t) ~-A(s)f g (s) + B(s)s df ^ t) 



dt --^'-y^' ■ v-/- ds 
For the case of strongly decreasing PTF the exponent under the integrals for the functions 
A(s) and B(s) can be expanded as 



A(s) 
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dqq 2 W(q); B(s) 



dqq 2 W'(q) 



(15) 
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Then the simplified kinetic equation for the case of short-range on q- variable PTF (non- 
equilibrium, in general case) reads 



df(s,t) 



A s 2 f(s) + Bqs 



df(s) 



dt ds 

where A = -1/6 / dqq 2 W(q) and B = -1/3 / dqq 2 W'(q). 
The stationary solution of Eq. ( fT4l) reads 

, A(s>) 



(16) 



f(s, t) = Cexp 



ds 



s'Bis' 



Cexp 



A s 2 
' 2B n 



(17) 



The corresponding normalized stationary momentum distribution is given by 

B p\ 



f(p) 



n 9 bT 



(2ttA)) 3 / 2 6XP L 2A 

Therefore, in Eq. ( ITT]) the constant C = N, where N is the total number of particles in 
the system, which undergo diffusive motion. Equation ( !T6|) and this distribution are the 
generalization of the Fokker-Planck case for normal diffusion in a non-equilibrium situation 
with strongly decreasing kernels W(q), W'(q) , when the prescribed PTF function ^(q, p) is 
determined, e.g., by some non-Maxwellian distribution of the small particles /&. To show this 
in an alternative way let us take the Fourier transformation of ffTTT) and the corresponding 
coefficients A and B n : 



(18) 



df(p,t) 



, d 2 /( P ,t) 0(p«/(M) 

-A.Q 7T- o Bq- 



dt u dp 2 dp a 

We then arrive at a Fokker-Planck type equation with friction coefficient (3 



(19) 



— B and 

diffusion coefficient D = —A /M 2 . In general these coefficients (Eq. f|T5|) ) do not satisfy the 
Einstein relation. 

In the case of equilibrium ^-function (e.g., /& = /°, see above) the equality W'(q) = 
W(q)/2MT b is satisfied. Then we find A(s)/sB(s) = A /B = MT b . In this case the 
Einstein relation between the diffusion and friction coefficients D = j3T/M is satisfied and 
the standard Fokker-Planck equation is valid. 

In the general case, however, the general equation ffTTj) . ( !T2|) and (JT3"]) have to be used. 
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III. DIFFUSION IN PLASMA SYSTEMS ON THE BASIS OF 
BOLTZMANN-TYPE COLLISIONS 

Let us calculate the PTF for the case of Coulomb collisions. The differential cross-section 
for the Coulomb scattering da /do equals 

da(x,u) f 2 . 2 x2 1 



(Z^n^) 2 - Ax = 2 (20) 



do sin' 2 , ; 

where \ = arccos(l — q 2 /2fi 2 u 2 ). Then 

^(q,P) = ^ r duu (Ze 2 /2^) 2l -^f b (u 2 + v 2 - q - V /^) 

8Z 2 e 4 7r f°° , t . 2 2 
= — / duuf b {u 2 + v 2 -q- v/fj). (21) 

It is necessary to stress that in the case of Coulomb interaction the general equations 
(llip . (114l) are applicable not only for diffusion of heavy particles in a light particle medium, 
but for arbitrary mass relations. The reason for this statement is the typical small transfer 
of momenta in the Coulomb systems. 

Let us calculate now the coefficients A a (p) and B a p(p) to compare the results with the 
linearized Landau kinetic equation, in which these coefficients depend on p. This implies 
that for the Coulomb interaction the expansion by (q • p) has to be performed at finite p. 

At first we consider the approximation in the spirit of the usual Fokker-Planck approach. 
Eqs. dHJ),© lead to the expressions 

MP) = J ^W^(q,p) = y J d*qq 2 w B c ° ul (q), (22) 
B a p{v) = \J ^qq^w^i^p) ^ l -J d 3 qq a qpwC° ul (q), (23) 



where 



8 Z 2 e 4 7r 

<? 5 Jq/2fi 



™ C B° u \q) = I duuf b (u% (24) 



w B Coul (q) = duufi(u% w B c ™\q) = --^ I duuf b {u 2 ) (25) 

and f b (y) = dfb(y)/dy. The procedure which we used here implies that the long tails of the 
functions w B oul (q) and w' B Coul (q) are absent. It is easy to see that the expressions (124|) . (|25|) 
in the limit of small q (the lower limit of the integrals in these equations is taken equal to 
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zero, which corresponds to the Landau small-g expansion) are appropriate in the Fokker- 
Planck equation to the Landau approach for the kinetic equation for plasma. In this case 
the coefficients A a (p) and B a p(p) read 

A I \ ~ P a f j3 2 Caul I \ no ]-G'K 2 Z 2 e' i p a ( q m ax\ t / oc x 

A a (p) = -jd qq w B (g) = In ( — j J, (26) 

B<# = *-fj d^w c B -\ q) - In (^£) J, (27) 

J = r duu f b (u 2 ) = n ^ (28) 

Jo (2^)3/2^ V ' 

Therefore, one can rewrite A a = —p a Vi ei where v ie is the characteristic frequency friction 
ions on electrons: 

4 v ^FZ 2 e 4 n/i 1 / 2 ^ (q max \ 



■limn 



3MT 3/ 2 y qn 

The corresponding friction force per unit volume F ie is equal Fi e = n,iM\Jv ie , where Fi e is the 
relative velocity of the electrons and ions [27]. In fact, the divergence at large q handled by 
a cut-off are not an artefact. This becomes clear when calculating the equilibrium function 
lOg " ' °(q, p) more accurately, without expansion on small values of q. 

For the equilibrium distribution function f®(u) = n e (fi/27iT) 3 / 2 exp(-fiu 2 /2T) the PTF 
function reads 



w 



Caul, 0/ 
B \ 



An Z 2 p 4 7r f°° 
q, p) = J^fJJL ex V [-^v 2 - q ■ v//*)/2T] / du 2 (fi/2nTf 2 exp[-/i M 2 /2T] 



q/2/i 

exp[-n(v 2 - q • y/(i + q W)/2T](30) 



An e Z e fj, / | . -2 .,./,. ..2/1.. 2 



^•°(q,p) = Jg^ g5 exp[-M, 2 - q ■ v//x + 9 2 /V)/2T], (31) 

Using the Fokker-Planck approximation for the coefficients A a (p) and B a p(p), and 
Eqs. we find 

MP) - | / 0) = gfgg£ * , (32) 

B«e = S fJ ^S"* "(a 0) = g^gg** /,, (33) 
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where 

J x = f d 3 qq- 3 exp[-q 2 /8/j,T] = 4vr / — exp[-q 2 /8{iT] ~ Air [ — exp[-q 2 /8{iT] 
J Jo Q ' ' Jg min Q 

= 2ir ^exp[-C] = -2-nEi{-q 2 /2^T)\Z m * -2irEi{-q 2 mm /2^T).{U) 
J sj)r s 

We can suppose that the minimal momentum transfer q m in is determined from the equality 
q 2 nin /2fiT = r m i n /r max . According to the Landau theory for a weakly interacting plasma 
r m in/r m ax = Ze 2 /Tr D < 1 for Ze 2 /hv T > 1, or r min /r max = h 2 /2fiTr D < 1 for the 
opposite inequality Ze 2 /Tivt <C 1. Here ro is the Debye radius and v? = y/T//j, is of the 
order of the thermal velocity In our approach the cut-off for the small momenta is satisfied 
automatically and corresponds to the second inequality (the "quantum" case). 

For a weakly non-ideal plasma this means a cut-off at the minimal momentum q m in = 
h/rn- Then 

Ji = -27rC + 47rln (— j ~4vrln (— — j , (35) 



where r^ in = fi 2 /2jxT and C ~ 0.577 is the Euler constant. 

It is easy to verify that W'{q) = W{q)/2MT and W"{q) = W{q)/AM 2 T 2 (in the case 
under consideration W(q, p) = Wg° u1, °(q, p)). Therefore, for the equilibrium case the usual 
Fokker-Planck equation for heavy particles (ions or dusty particles in dusty plasmas) is, 
naturally, valid with a good accuracy, owing to the exponential convergence of the integrals 
in the coefficients A(s) and B a (s) is provided at high values of q. The term with W" in 
Eq. ((Tj) is negligible, according to the above general statement. However, for small q the 
coefficients A and B a have the logarithmical divergence typical for Coulomb interaction 
because W ~ 1/g 5 just like W'(q). As follows from Eq. (I2TI) . this divergence not only exists 
for equilibrium, but for an arbitrary distribution function The simplest physical way 
to avoid this divergence is to cut the integrals for A and B a in Eqs. (TT2]) . (fl3|) for small q 
by the Debye radius 1/ro, following the well known the Landau procedure. We are more 
interested to find examples for non-exponential behavior of W, which may occur, e.g., for 
some specific non-Maxwellian distributions ft,. 
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IV. MODELS OF ANOMALOUS DIFFUSION FOR COULOMB INTERACTION 

Now we can calculate the coefficients for models of anomalous diffusion in plasma-like 
systems. 

At first we calculate the model of a Coulomb system with two species of particles with 
masses /i and M ^> /i. Let us suppose that in the model under consideration the small 
particles are described by a prescribed stationary distribution f b = n b (j) b /ul (where (f) b is 
the non-dimensional distribution, uq is the characteristic velocity for the distribution of the 
small particles) and £ = {u 2 + v 2 — q • v//j,)/uq. 

W«(q, P) = ^ di ■ MO- (36) 

U 0Q V(q2/4 jU 2 + „2_ q . v / /Lt )/ u 2 

First, let us consider a power type distribution 4>b(0 — C/£ 7 (7 > 1) 



^q.p = 5 — 71 d& = 5 7> 37 

u g 5 (1 - 7) so « g 5 7-1 



where £ = (? 2 /4yU 2 + u 2 — q ■ v//x)/wq. 

For the case p = the value £0 — ^ £0 = <? 2 /4/i 2 Uo and we arrive at the following expression 
for anomalous = W^: 

87rZ 2 e 4 n fe C(g/2^ ) 2 -^ 2 2 ^ 1 7rZ 2 e 4 n fe ^- 3 / i 2 ^ 2 C 
^ a(q ' P " 0) " — ^ 7^1^ " (7 - l)g*H-3 • ( 38 ) 

To determine the structure of the transport process and the kinetic equation in velocity 
space we have to determine also the functions W'(q) and W"(q). 

If p 7^ we have to use the full expression £ = (g 2 /4/i 2 + p 2 /M 2 — q • p/M/^/ii 2 , and it 
derivatives on q • p at p = 0: £q = —1/MfjLU^ and £q = 0. Then 

W q, p = -T7 — W q, p = +1 . 39 

For p = (£0 — >• £0) we obtain the functions 

o2 7 +3 2 7 -l 2 7 -3 7 2 4 n o2 7 +5 2 7 2 7 -3 72 4 

W ~ Mq*r+ B ' M) ~ M 2 q 2 ^+ 7 ' { } 

For the function A(s), according to Eq. (fl2|) . we find 

sin(qs) 



A(s) = Air 



POO 

/ dqq 2 


sin (qs) 


Jo 


qs 



W{q)=4*C a J~ 'dg^Lj 



- 1 



(41) 
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Comparing the reduced equation (see below) in velocity space with the diffusion in co- 
ordinate space (27 + 3 <H- a and W(q) = C a /q 2l+3 ) we can conclude the integral in the 
right-hand side of Eq. ( 14"TT) (3d case) converges if3<27 + 3<5or0<7< 1. The 
inequality 7 < 1 implies convergence for small q (q — > 0) and the inequality 7 > implies 
convergence for q — > 00. Likewise for the integral in B(s): 

sin(qs) 



4tt , . 2 



B ( s ) = -5- / d M 



s 2 



cos (qs) — 



W'(q); (42) 



qs 

convergence of B(s) exists for small q if 7 < and for large q — ?■ 00 if 7 > — 3/2. Again, it 
is easy to show that the term with W" can be omitted. 

Therefore, for convergence of A and B for a large g we require convergence for A, which 
implies 7 > 0. For convergence for small q it is sufficient to have convergence for B, implying 
7 < 0. Therefore, for the case of pure power behavior of the function /&(£) convergence is 
absent. It is also clear that the function /&(£) = C/tp (7 > 1) cannot be normalized. 
However, for anomalous diffusion in momentum space in reality the convergence for small q 
is always obtained, e.g., by a finite value of v or by a change of the small g-behavior of W(q) 
by screening (compare with the examples of anomalous diffusion in coordinate space [13]). 
Therefore, in the model under consideration, the "anomalous diffusion in velocity space" 
for a power behavior f b {u) (and as a consequence with a power dependence of W (q) and 
W(q)) on large q exists if for large q the asymptotic behavior of W(q — > 00) ~ l/q 2j+3 with 
7 > 0. At the same time the expansion of the exponential function in Eqs. (Tl2l) . (TIBl under 
the integrals, leading to the Fokker-Planck type kinetic equation, is invalid for power-type 
kernels ^(q, p). 

As an example of the above statements, let us consider the Cauchy-Lorentz-like distribu- 
tion for the function f b (r = 3): 



_ „ ^0 

'n 2 (u 2 + Vq) 

Then we find 



f b {u 2 ) = n b ^ /n 2 ; (43) 



Wa(q, p) = 5 — / w / fe (w 2 + v - q • v/fi) - 

Q Jq/2n 

Av 3 Z 2 

1 • e 

where £0 = (g 2 /4/i 2 + v 2 — q • v///)/uq and 



/- 00 ,^ 1 4^Z 2 e 4 rvr , I 

x / ^ = i 1 o ~ arctan£ \ , (44) 



s Av Z 2 e 4 1 . . 

W ^ = WWTJ- (45) 
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For large q the functions flUD and (l4"5j) tend to W (q, p) ~ 16fQyU 2 Z 2 e 4 /7rg 7 and W' a {(\, p) ~ 
64fQ/i 3 Z 2 e 4 /M7rg 9 . For small q convergence of the coefficients A(s) and B(s) cannot be 
obtained since these functions are determined by the expressions W a (q, 0) and W^(q, 0). 
However, this problem can be avoided by using a cut-off of the respective integrals (JUJ) 
and (j32J) at small q or by modification of the distribution (j4"5|) at small q (in the spirit 
of the respective cut-off for anomalous diffusion in coordinate space [12]). For large q the 
Cauchy-Lorentz-type of distributions have long tails, thus leading to anomalous diffusion. 

Let us now consider the formal general model for which we will not connect the functions 
W(q) and W'(q) with a concrete form of W(q,p). Therefore we consider the problem 
suggesting some behavior of the function W(q,p), but not on the level of the distribution 
function In general the functional VT(q, p) is unknown. In this case one can suggest that, 
independently one from another, the functions W(q), W'(q) and W"(q) possess a power- type 
g-dependence for a large q. 

As an example, this dependence can be taken as the power type for two functions W(q) = 
a/q a and W'(q) = b/q 13 , where a > and /3 are independent. Then, as follows from the 
consideration above, convergence of the function W exists if 5 > a > 3 (for asymptotically 
small and large q, respectively). For the function W'(q) the convergence condition is 5 > 
f3 > 2 for asymptotically small and large q, respectively. 

Finally for the function W"(q) the convergence condition is 7 > r] > 5 (for asymptotically 
small and large q, respectively). In this case the terms with W" can be omitted (for the 
same reasons as above). 

For this example the kinetic equation Eq. pip reads 

= P 0S -3 /(s , t) + /- 5 P lS ,A/(s, t), (46) 

where 

sin C 



P = 4na / d(( 



-2-a 



- 1 



Pi = Aub / dCC 



oo 

2-/3 



cos C — 



(47) 



(48) 



C 

Taking into account the isotropy in s-space we can rewrite Eq. fj46|) in the form 

= PoS ^f(s,t) + s^P^f^t), (49) 
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Naturally, Eqs. ( H6|) and (jUJ) can be formally rewritten in momentum (or in velocity) 
space via the fractional derivatives of various orders (see below). Therefore, as is easy to 
see, for the purely power behavior of the functions W(q) and W'(q) the solution with the 
convergent coefficients exists for powers in the intervals mentioned above. We establish that 
the universal type of anomalous diffusion in velocity space in the case under consideration 
exists if 5 > a > 3, 5 > > 2. This universality seems similar to the universality of the Levy 
distribution in coordinate space, where the power a of the dependence of the PT-function in 
coordinate space, W(q) ~ C/q a on the displacement q lies in the interval < = a — r < 2 
(r is the dimension of the coordinate space). 

As is easy to see the stationary solution of Eqs. (fl9l reads 



f(s) = C'exp 



(50) 



(a-/3 + 2)P 1 _ 

Of course, the general description above is also valid for the more complicated functions 
W and W, possessing a non-power dependence on q at small q and an asymptotical power 
dependence on q at large q. In this case the limitations for convergence are connected only 
with large values of q, namely it is enough to provide the inequalities a > 3 and (3 > 2. 
Simple examples of such type of PT-functions are (in analogy with anomalous diffusion in 
coordinate space [13]): 

W = ^(q)= 1 "y ) . (51) 

The correspondig kinetic equations in these cases cannot be written in partial derivatives 
and evolution of the system has to be described by Eq. (jTT|) . or for the isotropic case by 
( JT4"j) . If external forces are present, they have be included in the usual way in the left 
side of Eq. (jTT]) . Physically, this type of PT-function behavior can appear, in particular, 
for the case of a turbulent plasma, when the development of some instability can create 
a strong chaotic electrical field or irregular chaotic motion of one sort of particles with a 
prescribed non-Maxvellian distribution function. In such a turbulent plasma scattering with 
large transferring momenta can play a crucial role. 
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V. REPRESENTATION IN MOMENTUM SPACE AND CONNECTION WITH 
THE FRACTIONAL DIFFERENTIATION APPROACH 



As mentioned, in general Eq. ( II lj) cannot be written in terms of fractional differentiation. 
It confirms that the approach to anomalous diffusion based on Fourier-transformation of the 
PT- functions, in the form applied in this paper (see also [13], [19]) is a more general way for 
the problems under consideration. 

However, for the purely power-type dependence of the functions W(q) = a/q a and 
W\q) = b/qP, where a, (3 are independent and satisfy the inequalities 5 > a > 3 
and 5 > (3 > 2, Eq. ( 1461) is appropriate and can be represented after inverse Fourier- 
transformation in the following form (with the fractional derivatives) 



dt 



P AV(P, t) + P\A X (x + Ap^) /( Pj t ), (52) 



where v = [a — 3) (2 > v > 0); A = — 5 (2 > A > 0). Here we introduced the fractional 
differentiation operator in the momentum space A A /(p,t) = J dss 2X exp(—ips)f(s,t). 

Let us consider now formally a specific particular model of anomalous diffusion, for which 
we assume a structure of the PTF H^q, p) with a rapid (say, exponential) decrease of the 
function W'(q). Therefore, the exponential function under the integrals in the coefficients 
B(s) can be expanded, implying B(s) = B (or /3 = 5 and B = Pi in the notations of 
Eq. (1461)). At the same time the function W(q) = a/q a has a purely power dependence on 
q. The kinetic equation Eq. (H91) then reads 



rf/(S ' t] - P s a - 3 f(s, t) + B oSl ^-f(s, t), (53) 



at as 
or in momentum space, according to ([52 



= P A v f(p,t) - Bo-^foffat)]. (54) 

Eq. f )54|) is similar to the corresponding equation in [22], where a model of the Langevin 
equation with a constant friction frequency u = —B has been considered. In [23] a similar 
model with v = 3/2 has been applied to estimate the fusion rate in a hot rarified plasma. 

The stationary solution of Eq. ( 1531) is: 



f(s) = C'exp l-^f]. (55) 
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The corresponding distribution in p-space is proportional to the Levy-type distribution 

W(y,a') = (2y/n) / °° dtt sin(yt) exp(-t a '): 

f(p) — C' I d 3 sexp(— ips)exp [— — ^ ] = — — / dsssin{ps)exp [— ], (56) 
J ' vB§ p J vB 

with y ee p (vBq/Pq) 1 / v and a' = za As an example for the case v — 1 we find /(p) 

= Eo[(p 2 + 4P Vi? 2 )]2- (57) 

In the case v — 1 the long tail of the distribution is proportional to p~ 4 and the distribution 
f{p) corresponds with the Cauchy-Lorentz distribution. Normalization of the distribution 
f(p) leads to the value C = n/(27r) 3 , where n = N/V is the average density of particles 
undergoing diffusion in velocity space. A similar appoaach can be taken for other types of 
anomalous diffusion in velocity space. 



VI. CONCLUSIONS 



In the present paper the problem of anomalous diffusion for plasma-like systems in mo- 
mentum (velocity) space is consequently analyzed. A new kinetic equation for anomalous 
diffusion in velocity space has been derived recently in [19], without suggesting any sta- 
tionary equilibrium distribution function. We applied this equation to a system of charged 
particles with different masses to describe diffusion of heavy particles (ions, charged grains) 
in the surrounding light particles (electrons for the electron-ion plasma, electrons and ions 
for dusty plasmas). The distribution of the light particles can be non-Maxwellian, which is 
the cause for appearance of the long tails in the probability transition function. Conditions 
of convergence for the coefficients of the kinetic equation have been derived for a number 
of particular cases. It is found that a wide variety of anomalous processes in velocity space 
exists. 

In general the Einstein relation for such a situation is not applicable, because the sta- 
tionary state may be far from equilibrium. For the case of normal diffusion the friction 
and diffusion coefficients have been found explicitly for the non-equilibrium case. For the 
equilibrium case the usual Fokker-Planck equation in plasma is reproduced as a particular 

case. 
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